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ABSTRACT 



It is well known that the mass function for halos in CDM cosmology is a 
relatively steep power law (n(M) oc M-2) for low masses, possibly too steep to be 
consistent with observations. But how steep is the galaxy mass function? We have 
analyzed the stellar and gas mass functions of the first massive luminous objects 
formed in a ACDM universe, as calculated in the numerical simulation described 
in Gnedin (2000ab). Our analysis indicates that, as suspected, the relationship 
between luminous or gaseous matter and total (primarily dark matter) mass is 
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not a linear scaling relation — i.e., light does not trace mass. While the dark 
matter mass function is steep, the stellar and gas mass functions are flatter for low 
mass objects. The stellar mass function is consistently flat at the low mass end. 
Moreover, while the gas mass function follows the dark matter mass function until 
reionization at 2; ~ 7, between z = 1 and z = 4, the gas mass function also flattens 
considerably at the low mass end. At ;2 = 4, the gas and stellar mass functions are 
fit by a Schechter function with ctRi— 1.2±0.1, significantly shallower than the 
dark matter halo mass function and consistent with some recent observations. 
The baryonic mass functions are shallower because (a) the dark matter halo 
mass (differential) function is consistent with the Press-Schechter formulation 
at low masses n{M) oc and (b) heating/cooling and ionization processes 

appear to cause baryons to collect in halos with the relationship Mf, oc at 
low masses. Combining (a) and (b) gives n{Mij) oc Mj^ , comparable to the 
simulation results. Thus, the well known observational fact that low mass galaxies 
are underabundant as compared to expectations from numerical dark matter 
simulations or Press-Schechter modeling of CDM universes emerges naturally 
from these results, implying that perhaps no "new physics" beyond the standard 
model is needed. 

Subject headings: cosmology — galaxies: abundances, dwarf, formation 

1. Introduction 

Recent interest in modifying the standard ACDM theory of the origin of structure in the 
universe has been motivated largely by observations on sub-galactic scales. In particular, the 
most straight-forward interpretation of the results of theoretical calculations appear to be 
in conflict with the observed characteristics of spheroidal galaxies (i.e., bulges of spirals and 
small to moderate ellipticals) with respect to their individual physical properties, luminosity 
function, and clustering properties (Sellwood and Kosowsky 2000, Klypin et al. 2000, Navarro 
and Steinmetz 2000). These facts have led some authors to turn to alternative models 
such as quantum isocurvature (Peebles 1999ab), Self- Interacting Dark Matter (Spergel and 
Steinhardt 2000), and Warm Dark Matter (Bode et al. 2001). 

However, while measurements such as those of the cosmic microwave background ra- 
diation can be directly computed from theory in a relatively straight-forward manner, the 
abundance and clustering of luminous galaxies cannot at the present time be directly calcu- 
lated with any certainty. Such calculations are hampered by the classic problem of relating 
light to mass. The simplest way to relate dark matter density to galaxy counts is to as- 
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sume a proportionality in the fluctuations. This has led to the invention of the "spatial bias 
parameter 6." With the increasing sophistication of both theory and observations, we have 
discovered that bias must depend on both length scale and galactic type (because different 
types of galaxies have different clustering properties) and possibly on both local density and 
temperature (Blanton et al. 2000a). 

But what is the relation between dark matter halos and galaxies as a function of mass? 
Is it plausible to assume that each dark matter halo contains a galaxy having stellar mass 
and light proportional to the dark matter mass? The earliest phenomenological treatments 
explored this possibility, but obscrvationally it is clearly false. For example, the most massive 
(~ IO^^Mq) identifiable dark matter halos — found from large angle arcs and multiply 
imaged galaxies — correspond to clusters of distinct galaxies and do not contain supergiant 
galaxies of corresponding baryonic mass (10^^ — IO^^Mq). 

For dark matter halos less than or comparable to ~ IO^^Mq, the prevailing belief is that 
they will typically contain one "galaxy" and several satellites, but whether this description is 
accurate (and informative rather than definitional) is not known very securely. In the semi- 
analytic treatment of galaxy formation (Somerville and Primack 1999, Lobo and Guiderdoni 
1999, Chiu and Ostriker 2000) halos are identified with galaxies by a set of prescriptions 
that, overall, appears to correspond to reality. But, hydrodynamic simulations tend in this 
mass range to find several galaxies orbiting within a common halo (Cen and Ostriker 2000). 

For still lower mass galaxies (and presumably lower mass halos), we know that the 
mass/light ratio becomes large, implying that galaxy formation was inefficient for small 
systems. Thus, we are led to ask, could the physical bias against low mass systems naturally 
produce a mass function for galaxies significantly flatter than that of the dark matter halo 
mass function? If so, some of the apparent discrepancies between CDM cosmologies and 
observations might disappear. 

It may be that some of the other discrepancies, such as those that depend on lensing (Li 
and Ostriker 2000), shape (Moore et al. 1999b, Miralda-Escudc 2000, Navarro and Steinmetz 
2000), or small scale potential fluctuations (Moore et al. 1999a), will remain and require 
serious reexamination of the CDM paradigm. But the driving reason for reexamination — 
the much smaller than expected number density of small galaxies — might simply be a result 
of well understood physical properties. This has been suggested already (Silk 1977, Lobo 
and Guiderdoni 1999) on the basis of semianalytical arguments. The purpose of this paper 
is to see if the mass function of low mass galaxies ("faint" galaxies) is significantly fiatter 
than that for low mass halos, and the evidence from the hydrodynamic simulation examined 
(Gnedin 2000ab) is yes. 
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The organization of this paper is as follows. In § 2 we describe the hydrodynamical 
simulation used. In § 3 we show our analysis of the mass functions of the bound systems in 
the simulation. We discuss the why the galaxy mass function is relatively flat in § 4, and 
conclude in 5 5. 



2. The Hydrodynamical Simulation 

The simulation code used here is based on Smoothed Lagrangian Hydrodynamics and 
is fully described in Gnedin (2000a). It is an extension of that reported on in Gncdin (1995, 
1996), Gnedin and Bertschinger (1996), and Gnedin and Ostriker (1997). In brief, this 
simulation contains 128^ dark matter particles, a baryonic mesh of the same size, has a 
spatial resolution of 1.0 kpc//i, and a box size of 4^ {Mpc/h)^. The baryonic mass resolution 
is IO^'^Mq. In addition, the code includes heuristic star formation, with the creation of 
stellar "particles" whose dynamics only depends on gravity. Finally, reionization by these 
"stars" is modeled with the local optical depth approach which is able to approximately 
follow three-dimensional radiative transfer. Because the present day nonlinear scale is larger 
than the box size, the simulation can not be continued to the present, and was stopped at 
a redshift z — 4. Bound objects at each epoch are identified through use of the DENMAX 
(Bertschinger and Gelb 1991) algorithm. The simulation run analyzed here is the same as the 
"production run" analyzed in Gnedin (2000ab). The cosmological parameters are Qq = 0.3, 
Qb = 0.04, Oa = 0.7, h = 0.7, n = 1, and as = 0.9 (COBE normalized). 

Using results of a similar code, Gnedin and Ostriker (1997) showed that ionization is 
largely complete for a standard ACDM model at z = 7. In addition, Gnedin, Norman, and 
Ostriker (2000) recently showed that when the effects of gas cooling and star formation are 
included, massive bulges are able to collapse to densities exceeding that of the dark matter, 
resulting in structures that resemble observed objects in their sizes, shapes, and density 
profiles. 

3. Mass Functions of Spheroidal GalzLxies 

Here, we examine the mass functions of the dark and baryonic components of the bound 
objects formed in our hydrodynamical simulation. From the simulation results, we calculate 
the mass function ni{> M,), the number of objects with mass in component i > Mj, where 
i =dark matter (d), gas (g), or stellar particles (s). If dark matter, gas, and stellar particles 
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all trace each other, then 

ni(> Mi) ^ nj(> Mj) (1) 

with the mass scaling 

Mr _ h ... 

where fi is the fraction of the mass of the universe in species i. The dark matter fraction 
fd = Qd/^o = 0.867 is universal. Baryons are either in the form of gas or stellar particles, so 
that fg + fs = 0.133. The fraction of mass in stellar particles is given in Table 1 as a function 
of redshift. At all times, less than 4% of the baryons are in the form of stellar particles, so 
fg remains approximately constant. 

The results of the comparison of the mass functions are given in Figures 1-3. The dark 
lines are measured from simulations, whereas the gray lines are various fits. The thin grey 
lines are the dark matter mass functions scaled via equations (l)-(2). The thick gray lines 
are fits with a Schechter function $(M) oc m" exp{—m/ m^)dm for the stellar {z — 7 and 
z — A) and gas {z — A only) components. The fits consistently give values of o; — 1.2 ± 0.1 
when fit for objects with M, > IO^-^Mq {z < 7) or Mg > 10^-^Mq ~ 4). These mass 
cutoffs were chosen due to numerical resolution considerations (Ix and 10 x baryonic mass 
resolution for stellar mass and gas mass respectively). The exponential cutoff may not 
be very accurate, since there is a lack of large scale power due to the small box size of the 
simulation. 

As an additional check, we compared the bound objects' dark and stellar masses as 
shown in Figure 4. Also shown are power laws Ms oc Mj and Ms oc M^. The motivation 
for this power law is the finding by Gnedin (2000b) that in low mass objects, heating 
consequent to reionization leads to a filtering of the baryonic mass resulting with such a 
relationship. The power law Mg oc Mj is a generally a good fit, while Mg oc clearly fails. 

Figure 4 also shows that the "tail" of the distribution around the Mj relationship is 
skewed so that for a given stellar mass Ms, the dispersion towards lower Md is greater than 
the dispersion towards greater M^. Even given this skewed distribution, there are several 
significant outliers in the tail. In particular, four bound objects with stellar mass > 10^Mq//i 
have significantly enhanced stellar mass fractions. For these objects, the fraction of mass in 
stellar particles is greater than even the universal baryon/dark matter mass ratio (shown in 
gray). 

Shown in Figure 5 is the gas mass as a function of the dark matter mass. While the 
stellar mass showed a clear ~ Mj dependence for all values of Mg for all relevant epochs, the 
gas mass (consistent with Gnedin 2000b) shows an evolving relationship. The gas starts out 
following the dark matter mass Mg oc M^. After reionization, a relation Mg oc Mj develops 
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at low masses, with the mass at which there is a transition to Mg oc increasing with 
time. 

Finally, shown in Figure 6 is the stellar mass as a function of the gas mass of bound 
objects. Because the gas and dark matter are closely correlated up until 2; ~ 7, at these 
epochs, the relationship between stellar mass Mg and gas mass Mg is similar to that between 
the stellar mass Ms and the dark mass M^. However, between z — 7 and z = 4, gas "ejection" 
from low mass halos is clearly evident. This reduction in the gas mass fraction was noted 
in the analysis of Gnedin (2000b). In some cases, the mass in gas is less than the mass in 
stellar particles — i.e., most of the gas in the halo has been converted to stars or ejected. 

4. Discussion: Why Is the Galaxy Mass Function Flatter Than the Dark 



The logarithmic slope a ~ — 1.2 of the galaxy mass function is significantly shallower 
than the dark matter halo mass function, particularly for low baryonic masses {Mg or Mg 
< 10^~^Mq//i). Therefore identification of low mass halos with low mass galaxies is incorrect. 
For instance, a one-to-one mapping of dark matter halos to galaxies normalized to high 
baryonic mass halos would result in an overestimate of the number of low baryonic mass 
systems by several orders of magnitude. In particular, this slope is comparable to that 
observed for low mass galaxies (Cross et al. 2000, Blanton et al. 2000b). 

The explanation of the value of a is easy to explain as a combination of the dark 
matter halo mass function and the relationship between baryonic mass and dark matter 
mass in halos. The dark matter mass function follows the steep slope expected from the 
Press-Schechter (Press and Schechter 1974) approximation 



where n{Md) is the comoving number density per unit mass, po is the mean mass density 
today, 5c is the critical linear density for collapse, a is the rms mass fiuctuation, and we are 
considering small scale fiuctuations where 5c cr and cr is a very weak function of mass. 
From Figures 4 and 5, is it clear that 



Matter Mass Function? 
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d 1 



(3) 



M, oc M^^ 



(4) 



for z <1 and 



Mg OC Mt 



(5) 
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for z ~ 4 at low masses. Combining equation (3) with equations (4)- (5) gives 

n{Ms) cx M^^/^ forz<7 (6) 
n{Mg) oc for^~4 (7) 

for low masses. Thus this simple analytic derivation explains the simulation results: the 
stellar mass function appears to be fiat hj z = 7 while the gas mass function evolves from 
simple scaling with the dark matter for redshifts of2; = 15~7toa gradual flattening to 
look like the stellar mass function from z ~ 7 — 4. 

Without additional analysis, it is difficult to determine whether the similarity in the 
two mass functions is (a) coincidental — two different mechanisms, (b) correlational — 
same mechanism of both, or (c) causal — stars — > gas. Clearly, though, the key to this 
question is understanding the origin of the ~ behavior for both stars and gas, and 
how they are related. Whatever its origin, the result reported above appears to be robust, 
since Nagamine et al. (2001) found the same flattening in the stellar mass function in a 
simulation using a completely different numerical treatment of gas (e.g., Eulerian rather 
than Lagrangian), radiation, and star formation. This work, using very different numerical 
techniques (but a comparable mass resolution) also found a ~ —1.2 for the baryonic mass 
function. 

Let us also consider the sensitivity to the star formation efficiency. If the star formation 
prescription underestimates star formation and all the remaining gas should have turned into 
stars hy z = 4, the stellar mass function would retain largely the same shape, but with a 
higher normalization. Only if the prescription overestimates star formation would the shape 
change significantly. In this case, the mass function would be even flatter, and ACDM could 
have the opposite problem of underproducing low mass galaxies! 

In fact, assuming that all the gas in halos at z = 4 turns into stars by 2; = (and 
neglecting the effects of merging) gives a "turnover" at a circular velocity of ~ 40 km/s. 
That is, the stellar mass function would follow the dark matter for Vc > 40 km/s, and 
would have a logarithmic slope of ~ — 1.2 for i>c < 40 km/s. This corresponds closely to the 
observed circular velocity distribution of galactic satellites as determined by Klypin et al. 
(2000). Thus, if additional gas ejection is signiflcant between z — A and z — 0, ACDM could 
indeed underproduce low mass galaxies! 

5. Conclusions 

We have analyzed the properties of the flrst bound luminous objects formed in a ACDM 
universe, as calculated in the numerical simulation described in Gnedin (2000ab). This 
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simulation is one of only a few full-scale, non-equilibrium hydrodynamical simulations to 
date which can individually resolve the components of spheroids in the low mass range 
comparable to low mass galaxies (see also Nagamine et al. 2001). 

We have found that the galaxy mass function at 2; = 4 is significantly flatter than 
the dark matter mass function, whether one considers "galaxies" identified by halos' stellar 

particles or gas mass. A flattened low mass end is found for the stellar mass function from 
z = 7 — A, and for the gas mass function at z ^ A. Taking into account the uncertainties due 
to the prescription for star formation, at 2; = 4 the stellar mass function in a ACDM universe 
is predicted to have a low mass logarithmic slope of a ~ —1.2 or flatter. The reason for these 
shallower baryonic mass functions is a combination of two effects: (a) the dark matter halo 
(differential) mass function is consistent with the Press- Schechter formulation at low masses 
of n{M) oc and (b) heating/cooling and ionization processes appear to cause baryons 
to collect in halos with the relationship Mi, oc Mj at low masses. Combining (a) and (b) 
gives n{Mb) oc comparable to the simulation results. Thus, the "underabundance" 

of low luminosity galaxies relative to halos emerges naturally from these results, implying 
that no "new physics" is necessarily required to understand the low mass end of the galactic 
mass function. However, because gravitational clustering is still evolving between z = 4, 
when the simulation ends, and the present epoch, it remains to be seen whether ACDM with 
the appropriate astrophysics can be reconciled with the observed clustering properties of low 
mass galaxies. 

Discussions with Ken Nagamine are gratefully acknowledged. This work was supported 
in part by NSF Grant AST-9318185. 
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Fig. 1. — Mass functions for bound objects at z = 15, before reionization. Shown are the 
dark matter (thin sohd), gas (thin dashed), and stellar particles (thick solid). Also shown 
(dashed grey) is the dark matter mass function with the mass scaled by the universal ratio 
of baryons to dark matter Qb/^^d- 
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Fig. 2. — Mass functions for bound objects at z = 7, immediately after reionization. In 
addition to the quantities in Figure 1, also are shown the dark matter mass function with 
the mass scaled by the ratio of stellar particles to dark matter ^g/^d (thin grey), and a 
Schechter function fit with a ~ —1.3 (thick grey) to the bound stellar systems. 
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M IM hL 



Fig. 3. — Mass functions for bound objects ai z — A. In addition to the quantities shown in 
Figure 2, also is shown a Schechter fit with a ~ —1.1 (thick dashed grey) to the gas mass 
function. 
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Md IM hL 

Fig. 4. — Scatter plot of stellar mass as a function of dark matter mass for bound objects at 
z = lb (open circles), z = 7 (small gray points), and z = A (small black points). The solid 
black line has Ms oc corresponding to the related finding by Gnedin (2000b) that M^, oc 
M^^^. The remaining three lines are linear relations where the constants of proportionality 
are (a) the universal baryon/dark matter ratio (gray) and (b) the global stellar mass/dark 
matter ratios z — 1 (short dash) , and z — A (long dash) . 
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Md IM hL 



Fig. 5. — Scatter plot of gas mass as a function of dark matter mass for bound objects 
at 2; = 7 (small gray points), and z — A: (small black points). The solid black line has 
Mg oc MJ. The linear relation (gray line) has a constant of proportionality equal to the 
global baryon/dark matter mass ratios. 
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Mg YM hL 

Fig. 6. — Scatter plot of stellar mass as a function of gas mass for bound objects at 2; = 15 
(open circles), z = 7 (small gray points), and z = A (small black points). The solid black 
line has M, oc M^. The remaining three lines are linear relations where the constants of 
proportionality are (a) unity (gray) and (b) the global stellar mass/gas mass ratios at 2; = 7 
(short dash), and z = A (long dash). 
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Table 1. Fraction of Mass in Stellar Particles 



Redshift 




15 


6 X 10-^ 


7 


9 X 10-^ 


6 


0.002 


5 


0.003 


4 


0.005 



